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Abstract 

We study the resummation of singular distributions for electroweak- 
boson production at large transverse momentum in hadronic collisions. 
We describe the factorization properties of the cross section near par- 
tonic threshold and construct the resummed cross section at next-to- 
leading logarithmic accuracy in moment space in terms of soft gluon 
anomalous dimensions. We present full analytical results for the ex- 
pansion of the resummed cross section up to next-to-next-to-leading 
order. Our results can be applied to W ^ Z , and virtual 7 production 
at hadron colliders. 
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1 Introduction 

The study of electroweak-boson production in hadron colliders provides use- 
ful tests of the standard model and estimates of backgrounds to new physics. 
For example, Wbb production is the principal background to the associated 
Higgs boson production, pp — >• bb)W, at the Tevatron [|l|. The com- 

plete calculation of the next-to-leading-order (NLO) cross section for W, Z, 
and virtual 7 production at large transverse momentum in hadron collisions 
has been presented in Refs. [0, ^ , following earlier results for the non-singlet 
cross section in Ref. [Q. 

The calculation of hard scattering cross sections near the edges of phase 
space (partonic threshold), such as electroweak-boson production at high 
transverse momentum, involves corrections from the emission of soft gluons 
from the partons in the process. At each order in perturbation theory one 
encounters large logarithms that arise from incomplete cancellations near 
partonic threshold between graphs with real emission and virtual graphs. 
These threshold corrections exponentiate and have been resummed explic- 
itly at next-to-leading logarithmic accuracy for a number of processes includ- 
ing heavy quark p, ^ 0, ||, H, |10|, dijet |]lT], |^, and direct photon produc- 
tion [|r^, n], using general techniques developed originally for Drell-Yan 



production |jT5|, |T6[. For a review see Ref. 



In this paper we study electroweak-boson production at large transverse 
momentum in hadronic collisions where these considerations are of relevance. 
Following Ref. |T^, we discuss in Section 2 the factorization properties of 
the single-particle inclusive cross section |T^, |, |TT], and identify its sin- 
gular behavior near threshold. We then proceed to resum the leading (LL) 
and next-to-leading (NLL) logarithms explicitly to all orders in perturba- 
tion theory. In Section 3 we provide full analytical results for the expansion 
of the resummed cross section at NLO and at next-to-next-to-leading order 
(NNLO). Our NLO expansion agrees near partonic threshold with the ex- 
act NLO calculations in Refs. 0, |^ while our NNLO results provide new 
predictions. 
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2 Factorization and resummation 

For the hadronic production of an electroweak boson V of mass my, where 
V ^W,Z, 7*, 

hA{PA) + hB{PB)^V{Q) + X, (2.1) 
we write the factorized single-particle cross section as 

X Eq ^f^I^A^^9}l^{s, t, u, Q, 1,^, (2.2) 

where Eq = (5°, is the parton distribution for parton / in hadron h, and 
a is the perturbative cross section. The initial-state coUinear singularities are 
factored into the 0's at factorization scale fxp, while /ir is the renormalization 
scale. 

At the parton level, the subprocesses for the production of an electroweak 
boson and a parton are 

qiPa) + giPb) q{Pc) + v{Q), 

qipa) + q{pb) 9ipc) + V{Q). (2.3) 

The hadronic and partonic kinematic invariants in the process are 

S={Pa + Pb?, T={Pa-Q)\ U={Pb-Q)\ S2 = S + T + U -Q\ 
s = [Pa+Phf, t^ipa- Qf, u^iph- Qf, S2^S + t + U-Q^, (2.4) 

where 5*2 and S2 are the invariant masses of the system recoiling against the 
electroweak boson at the hadron and parton levels, respectively. S2 = {Pa + 
Pb — QY parametrizes the inelasticity of the parton scattering: for one-parton 
production, S2 = 0. Since Xi is the initial-parton momentum fraction, defined 
by Pa = XiPa and Pb = X2PB, hadronic and partonic kinematic invariants 
are related by s = X1X2S, t — — Xi{T — Q^), and u — Q'^ — X2{U — Q^). 

In general, a includes distributions with respect to S2 at nth order in ag 
of the type 

' In"'is2/Q') 

S2 



m 



<2n-l 



(2.5) 
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defined by their integral with any smooth function / by 



ds2 f{s2) 



S2 



Q\ \rris2/Q' 



ds 



S2 



■[f{s2)-fm. (2.6) 



These plus distributions are the remnants of cancellations between soft and 
virtual contributions to the cross section. 

We now consider the partonic cross section, with the colliding hadrons 
in Eq. ( p.2|) replaced by partons. To organize the plus distributions in a, 
we introduce a refactorization p, |Tl], [12| using new functions ■?/'/// and J 
that describe the dynamics of partons moving collinearly to the incoming 
and outgoing partons respectively, and functions H and S which describe 
respectively the dynamics of hard partons and of soft partons which are 
not collinear to ipf/f and J. This factorization is shown in Fig. 1 in cut 
diagram notation, showing contributions from the amplitude and its complex 
conjugate, with H = hh*. At partonic threshold we find the relation 



S2 
S 



u 



'l 



J S2 

X2)-^ 

S S 



= Wi 



-u 



+ W2 



-t 



+ Wj + Ws 



(2.7) 



where in the second line of the equation the weights w identify the contri- 
butions of the functions in the refactorized cross section. At fixed S2, the 
partonic cross section factorizes into 



Ec 



d^Q 



H / dwidw2dwjdwsipf^/f,{wi)ipf^/f^{w2) Jiwj) S{wsQ/ni 



S, 



-u 



Wj - Ws 



(2.8) 



If we take moments of the above equation, with the moment variable, we 
can write the partonic cross section as 



dS2 NS2/s do'fah- 

^^^^ 



s 



H I dw.e-'^^-^^j^/fSw,) 



dw2e 



-N2W2„ 



^h/hM I dwje-^'^'Jiwj) I diVse-'^^^^SiiVsQ/fiF) 



-Nws 



(2.9) 
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(a) (b) 

Fig. 1. (a) Factorization for V + jet production near partonic threshold, 
(b) The soft-gluon function S, in which the vertices c hnk eikonal hues rep- 
resenting the partons in the process. 



with A^i = N{—u/s) and N2 = N{—t/s) for either partonic subprocess. 

Taking moments of Eq. ( |2.2| ) with the incoming hadrons replaced by 
partons, and using the first line of Eq. (|2.7| ), we also have the relation 



X 



dS2 „_iVs2/sz7 ^^/aA-v(s2) 



Note that S2/S = S2/S up to quadratic terms in (1 — xi) and/or (1 — X2) 
and/or S2- Comparing Eqs. ( p.lOp and ( |2.9| ) we then may solve for the 
moments of the perturbative cross section a: 

E.'lM^Jj^^^ (2.11) 

The plus distributions in a produce, under moments, powers of In as high 
as In " iV. The LL corrections are included in ip/cj) and J, while the NLL 
corrections are included in '0/0) J, and S. 

The resummation of the A^-dependence of the jet and soft functions in Eq. 
( |2.11| ) depends on their renormalization properties |^ The factor 
is universal between electroweak and QCD-induced hard processes and its 



resummation was first done in the context of the Drell-Yan process [15, IG 
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It contributes an enhancement to the cross section, while the final-state jet 
J gives a negative contribution [jll], |12|. The UV divergences induced by 
factorization in the hard and soft functions cancel against each other since 
there are no additional UV divergences aside from those already removed 
through the usual renormalization procedure. The renormalization of the 
hard and soft functions can be written as 



(2.12) 



i=a,b 



where Zs is the renormalization constant of the soft function, and Zi is the 
renormalization constant of the ith incoming partonic field. Then S satisfies 
the renormalization group equation 



d 



d 



9 d 



/^T^ + S = -2(Rer5) 5, Tsig) = -^—Res^^^Zsig^e 



dg 



2dg 



(2.13) 

where is the soft anomalous dimension. We may then evolve S from the 
scale Q/N to fip as 



S l-^,a,{fil)] = S {l,a,{QyN')) exp j"'" ^ 2Rers (^.(/x'^) 



(2.14) 



The resummation of the A^-dependence of each of the functions in the refac- 
torized cross section, Eq. ( p.ll| ), in the MS factorization scheme, then gives 



Ec 



da 



fafb~*V 



(N) 



X 



i=a,b 

exp{E[^)(iV)} H(a,{^l 



IhMf^' )) -if^f,{Ni,as{fi' )) 



- / \ rQ/'^ da' / 

X 5(l,a,(QViV2)) exp Jf- 2ReT s [asifi 



,(2.15) 
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with = p'^/Q. The incoming parton-jet anomalous dimensions [|TT], |1^ 
defined through 



and 



/i 



d'iljf/f{N,Q/fi, 
d(j)f/f{N, fi'^,e) 



(2.16) 



(2.17) 



are given by 



7g 
75 



3„ as 
-Up — 

4 71 
4 TT ' 



7,9 = - IniV- 



IT 



793 



(2.18) 



for quark and gluon jets, respectively. /3o = (IIC^ — 2r2j)/3 is the one- 
loop coefficient of the (3 function, with Uf the number of quark ffavors. The 



exponents for the incoming jets are 0, |Tl], |12 



— I dz 



1 ^N,-l _ 1 



1 - Z 
1 



1 



(1-2)2 A 



A 



as {{2pi ■ O'A 



if:) 



a. 



((2p. ■ C)'(l - 



(2.19) 



At next-to-leading order accuracy in In(A^), we need u'^^^ = 2Cfas/'Tc and 
y4(j)(as) = Cf (as/iT + {K/2) (a^/vr)^^, with Cf = Cp (Ca) for an incoming 
quark (gluon), and K = Ca(67/18 - - 5n//9. 

The exponent for the final-state jet is |TT|, [T2[ 



1 ^N-l _ 1 



dz- 



1 - z 



(1--) dX 

(l-z)2 A 



A 



(J) 



+ B[j^ [a, ((1 - z)Q')] + B'(j^ [a, ((1 - zfQ')] } (2.20) 



where^ 



B' 



a. 



TT 



a. 



-^[-- Cp, B''=^Cp[H2u,)~l] 



TT 



^Note that Eq. (2.20) corrects the arguments of ag in B' in Eqs. (11) and (12) of 

Ref. in. 
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(a) 



TT 



' 4 



a. 



TT 



Ca [H2pg) - 1] 



(2.21) 



n] 



The i/j terms are gauge dependent. They are defined by z/j = {j3i 
where /?j = Pi\j2/ s are the particle velocities and n is the gauge vector, chosen 
so that Pi ■ C, = Pi ■ n for i = a,b [|l^. In our calculations we use Feynman 
rules for eikonal diagrams in axial gauge (resummation calculations have also 
been done in covariant gauge |T^ ). 

The soft anomalous dimensions are calculated explicitly by evaluating 
one-loop vertex corrections [T^. For the qg — > qV channel in the kine- 
matics CO) we find 1121, |T7| 



J- S 



+ C 



2 In (^^j -ln(4z/,„z/,J + 2 
In ( - \n{2ug) + l-TTi 



(2.22) 



The soft anomalous dimension for qq — > gV is [|T^, [17 

g.gj , 2 - 2m] + Ca 



J- S 



a 



Ztx 



In - \n{2vg) + l + Txi 

(2.23) 

Eqs. ( p.22| ) and ( |2.23| ) coincide with the corresponding soft anomalous dimen- 
sions for direct photon production [|12|, [T^ . Substituting Eqs. ( p.l9|) through 
( |2.23D in the resummed cross section, Eq. ( p^.lSj ), we see that at NLL accu- 
racy all the gauge- dependent terms cancel out in the exponent, for both the 
qg — > qV channel and the qq — > gV channel. 

We can rewrite the resummed cross section in a form which is more con- 
venient for the calculation of the fixed-order expansions. Using the renor- 
malization group behavior of the product HS from Eq. ( |2.12| ), 



^Jij-\n[H{^^)S{Q/{N^^))] 



and the relation 



J flu fi 



(2.24) 



(2.25) 
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with (3{as) = jJidln g / djj, = —f3oas/4:T[ + we can write Eq. ( p.l5| ) as 



d^Q 
X exp< ^ 



t=a,b L 



X exp{E[JN)} exp 
X S{l,as{QVN')) exp / 2ReTs {a^ifi 



(2.26) 



3 Next-to-next-to-leading order expansion of 
the resummed cross section 

In this section, we expand the NLL resummed cross section up to NNLO, in- 
vert back from moment space to momentum space, and perform a comparison 
with the NLO resuhs of Ref. @, ^. At the parton level, the subprocesses for 
the production of an electroweak boson and a parton are given in Eq. ( p.3|) . 

We can write the NLO corrections from Eq. ( p. 261) for qg — > qV in 
single-particle inclusive kinematics at NLL accuracy as 



Er 



da. 



NLO 
qg^qV 



d^Q 



a 



as(/i|) 



qg^qV 



IT 



{Cp + 2Ca) 



Hs2/Q') 



S2 



tu 





■ 1 ■ 




.S2. 



+ 5(S5 



Here the Born term is 



a, 



qg^qV 



aas{^l\)CF Mq ^l\T |2 , IP |2n 



s{N^c - 1) 



(3.1) 



(3.2) 



/ 
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\t s st 



with L and R the left- and right-handed couphngs of the electroweak boson 
to the quark hne, / the quark flavor and J2 / the sum over the flavors allowed 
by the CKM mixing and by the energy threshold. We choose for the L , R 
couplings the conventions of Ref. 0. The Born differential cross section is 
Eqda^g^qy I d?Q = af^^gv/ (5(s2). 

We have kept the factorization and renormalization scales separate. The 
expansion gives all [ln(s2/(5^) / ^2] + and [1/52]+ terms but only scale-dependent 
5(52) terms. 

Our one-loop expansion can be compared with the exact NLO cross sec- 
tion [0] in the proximity of the partonic threshold, S2 — > 0, i.e. with the sin- 
gular terms of the type ( |2.5|) in the contribution of the subprocess qg — > qVg 
to the cross section for qg — > qV. We find full agreement^ with the results 
of Ref. i. 

Then we expand the resummed cross section, Eq. ( |2.26D , for qg — > qV to 
two-loop order and compute the NNLO corrections in single-particle inclusive 
kinematics at NLL accuracy, including all the factorization and renormaliza- 
tion scale dependent terms. We obtain 



Er 



+ 



da, 



NNLO 
qg^qV 



3 



B 

qg^qV 



TT 



{Cp + 2CA)\{CF + CA)\n 




{Cf + 2Ca 



+ In 




I3q (Cf ^ 



{Cf + CaY In ( 1 + -^C'f + 2Ca{Cf + Ca) In , 




+ (C. + 2C.)(c.lnfeUc.ln g 



HS2/Q^) 



S2 



S2 



^Note that Gonsalves et al. use "A+" distributions, which are related to our "4- 
distributions by [ln"(s2/Q')/s2]+ = [^^'^ {s2 / Q"") / S2W + [1/ {n + I))\n^+\A/Q^)5{s2). 
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+ ln^ 



{Cf + Ca) 



In 



—u 



4 







■ 1 ■ 






.S2. 



^^{Cf + Ca) ln(^ 



») 


■ 1 " 




-S2- 



(3.3) 



Our expansion gives all \\.v?{s2/ Q"^) / S2\+ and [ln^(s2/Q^)/s2]+ terms but 
only scale-dependent terms for [ln(s2/Q^)/s2]+ and [l/s2] + . However we 
can derive all the [ln(s2/Q^)/s2]+ terms by matching with the exact NLO 
cross section^. Including the /x/j^iT'-dependent terms given above, the full 
[ln(s2/(5^)/'52]+ terms are 



cr 



qg^qV 



^^)' {^(C^F + 2Ca) [Br + Bl^uf + cr + Cfuf 



+B: 



qg {Lfafa - Rfafa) f{L f f - Rff) 

j:A\Lfu? + \Rfu\') 



+ 



{Cf + C a) In 



+ -Cf + Ca In 



tu 



+ {Cf + 2C^)9 In ^ + -K{Cf + 2Ca) - C2{Cf + 2Ca) 



c 



Ca 



tu 

7q^ 



ln(s2/g^) 



S2 



(3.4) 



with (2 = 71"^ /6, and with Br, Br, Br, Cr , and as given in the Ap- 
pendix of Ref. 0] but without the renormalization counterterms and using 
fA ^ ln(A/g2) = 0. 

Next, we consider the qq — ^ gV partonic subprocess, and compute the 
NLO corrections from Eq. ( p. 26 ) in single-particle inclusive kinematics at 
NLL accuracy. 



Er 



da 



NLO 
qq^gV 



d^Q 



a. 



B_ 

qq-^gV ' 



as{flR) 



TT 



{ACf - Ca) 
tu 



2CFln(!^]+{2CF-CA)\n 



S2 
00 

4 



1 

S2 



Eq. (2.12) of Ref. |Q] there is a typo: the term should be replaced by 
Sis2)Crnf + Cr- 
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In 



Cf In 



Here the Born term is 



a 



gg^gV 



7 + - + ^ 

t u tu 



tu 



In 



2, 



+ \R 



fbfa 



(3.5) 



(3.6) 



and the Born differential cross section is Eqdagg^gy/d^Q 



a 



gg^gV 



5(^2) 



The one-loop term can be compared with the singular terms of the type 
^ in the contribution of the subprocesses qq — > gVg and qq — > Vqq to 
the cross section for qq — > gV. We are in agreement with Ref. 0. 

The NNLO corrections from Eq. (|2.26|) for qq — > gV in single-particle 
inclusive kinematics at NLL accuracy, including all the factorization and 
renormalization scale dependent terms, are 



Er 



da 



NNLO 
gg^gV 



d^Q 
+ 



a. 



B 

gg^gV 



'a^(/i 



TT 



-(4CF-CAf 



S2 



^i5CF 



+ ln 



Cf ACf In 



2 



\+{2CF-CA)\n 



tu 



+ 4 {2Cf - Ca) In 



S2 

+ (4C^ - Ca 
tu 



In 



-I- In 



+ ln' 



1 2 



{^Cf - Ca) 



c. 



'/4 



Cf 


;-.,.( 






) 


(S)l 





ln(g2/g^) 

S2 

tu ' 

1 

S2 



S2 



3' 



3 + 



/3o 



1 

S2 



(3.7) 



Again, only scale-dependent terms for [ln(s2/Q^)/s2]+ and [1/^2]+ are in- 
cluded in Eq. (p^). The full [Iia{s2/Q^)/ 82]+ terms are found by matching 
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with the exact NLO cross section, 



a 



B 

qq^gV 



TT 



1 



{4Cf - Ca) \Bf + Cf + {Bf + 



Hi 



(Lfafa - Rfafa) ^ f f f " Rff) 



+ 



2Cpln 



4 



tu 



+ (2Cf - Ca) In 



+ 



+ (4Ci. - In ( ^ ) + (4C^ -Ca)^ + ^- U^Cf - Ca? 



In 



( tu 



ln(s2/Q^ 

S2 



(3.8) 



with Bf, Bf, Cf, and D^o) as given in the Appendix of Ref. § but 
without the renormahzation counterterms and using J'a = 0. 



4 Conclusion 

We have presented the resummed cross section at NLL accuracy for electroweak- 
boson production at large transverse momentum near partonic threshold. 
The expansion of the resummed cross section at NLO agrees with previous 
exact NLO results while the NNLO expansion provides new predictions at 
higher order. In future work we plan to study the numerical significance 
of resummation for W+ jet production. Related studies for other processes 
have shown reduced factorization scale-dependence, expected on theoretical 
grounds pOi, and increases over the NLO cross section. 
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